Can we ever distinguish between quintessence and a cosmological constant? 
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Many ambitious experiments have been proposed to constrain dark energy and detect its evolution. 
At present, observational constraints are consistent with a cosmological constant and there is no 
firm evidence for any evolution in the dark energy equation of state w. In this paper, we pose 
the following question: suppose that future dark energy surveys constrain w at low redshift to be 
consistent with —1 to a percent level accuracy, what are the implications for models of dynamical 
dark energy? We investigate this problem in a model-independent way by following quintessence 
field trajectories in 'energy' phase-space. Attractor dynamics in this phase-space leads to two classes 
of acceptable models: 1) models with flat potentials, i.e. an effective cosmological constant, and 
2) models with potentials that suddenly flatten with a characteristic kink. The prospect of further 
constraining the second class of models from distance measurements and fluctuation growth rates 
at low redshift (z < 3) seems poor. However, in some models of this second class, the dark energy 
makes a significant contribution to the total energy density at high redshift. Such models can be 
further constrained from observation of the cosmic microwave background anisotropics and from 
primordial nucleosynthesis. It is possible, therefore, to construct models in which the dark energy 
at high redshift causes observable effects, even if future dark energy surveys constrain w at low 
redshift to be consistent with —1 to high precision. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Ever since the discovery of cosmic acceleration P, Q , 
dynamical models of dark energy have been considered 
as serious alternatives to the cosmological constant (see 
for review). In the simplest dynamical models in- 
volving a scalar 'quintessence' field [J, |5|, 16|] rolling along 
a potential, the equation-of-state parameter w = p/p is 
time- dependent with w > —1. This contrasts with the 
cosmological constant, for which = —1 at all times. 

To test the nature of dark energy, a variety of astro- 
nomical techniques has been used to constrain w. We 
will not attempt to review the observational results in 
great detail here, but instead refer the reader to recent 
papers 0, 0, Q that describe results from combinations 
of type la supernova (SNIa) surveys, cosmic microwave 
background (CMB) anisotropics, galaxy redshift surveys 
and other cosmological data. These analyses are con- 
sistent with w = —1 to ^ 10% accuracy, depending on 
exactly which datasets and theoretical assumptions {e.g. 
curvature, evolution of w) are used. Although some data 
have suggested otherwise [lo|, [HI , there is a general con- 
sensus that the cosmological constant provides an excel- 
lent fit to the data and that there is no firm evidence for 
dynamical dark energy. 

Nevertheless, understanding the nature of dark energy 
is of such importance for fundamental physics that a large 
number of ambitious new surveys have been proposed to 
tighten the limit on w and its time evolution. These ex- 
periments include extensive ground and space-based sur- 
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veys of distant type la supernovae [ij, [13], the Planck 
CMB satellite [14], large galaxy surveys to measure accu- 
rately the baryon acoustic oscillation (BAO) [l5|,[l^ over 
a wide range of redshift, X-ray and Sunyaev-Zeldovich 
cluster surveys [13, [3, El] and ambitious ground and 
space-based weak lensing surveys 0, Hil. We will not 
discuss the relative merits of these methods here. Sum- 
maries of many of these projects can be found in [22l,[23|. 

Provided systematic errors can be kept under control, 
it may be possible to constrain w to an accuracy of a 
percent or so within the next decade. Either we will 
find definitive evidence of a departure from w = —1, in 
which case dark energy must be dynamical, or the obser- 
vational constraints will continue to tighten up around 
the cosmological constant value of w — —1. The problem 
that we wish to address in this paper is the following: if 
observations constrain w to be very nearly —1 to high 
accuracy over some range of redshift, does this imply a 
cosmological constant or can we always construct radi- 
cally different dynamical models that are consistent with 
the data? 

In this paper, we will approach this problem in a 
model-independent way, without making specific as- 
sumptions about the form of the quintessence potential. 
While there are many quintessence models in the litera- 
ture, some (loosely) inspired by string theory or super- 
gravity, we clearly do not yet have a fundamental theory 
of dynamical dark energy. Given our current poor un- 
derstanding of the physics of dark energy, quintessence 
models should be assessed primarily on their ability to 
explain observations. 

There have been many investigations of model- 
independent constraints on quintessence 0, [24|, |25|, [26|, 
[27| . However, our approach, based on 'energy' phase- 
space is quite different from these analyses. By inves- 
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tigating the attr actor dynamics in this phase-space, we 
show in an intuitively clear way why it may well be im- 
possible to ever distinguish between the cosmological con- 
stant and dynamically evolving models of quintessence. 



is analogous to the slow-roll parameter e in inflation. In 
the next Section, we explain how A can be used to gen- 
erate quintessence models stochastically. 



II. DYNAMICS IN THE ENERGY 
PHASE-SPACE 

A. Energy variables 

In this paper, quintessence potentials V{(j)) are gener 
ated using the phase-space approach introduced in 
and subsequently used by a number of authors [2^, 
to study the dynamics of speciflc quintessence potentials. 
To keep the discussion simple, we ignore any coupling of 
the quintessence fleld to matter, and we assume that the 
universe is spatially flat. In the phase-space approach, a 
quintessence model is represented by a trajectory in the 
plane of 'energy' variables (x, y) deflned by: 



X = 



7^' 



kVv 
7^' 



(1) 



where k = = vSttG and H = d/a is the Hubble 

parameter. These variables give the ratios of kinetic and 
potential energy of the quintessence fleld to the total en- 
ergy density. Energy conservation then requires: 



(2) 



where and are the matter and radiation density 
parameters. This constrains the quintessence trajectories 
to always lie within the unit circle in the energy phase- 
space. By restricting ourselves to V{<p) which are positive 
and monotonically decreasing, we only need to consider 
trajectories that lie within the flrst quadrant of the en- 
ergy phase-space. These assumptions will not affect our 
main conclusions and are discussed further in Section IVTl 
Using the Friedmann equations, one can express the 
trajectories as a function of redshift z: 
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dz 
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,(3) 



where the equation of state of the matter-radiation fluid, 
Wb, is during matter-dominated era, and 1/3 during 
radiation era. The 'roll' parameter A given by 



A = 



_}_dV 



(4) 



B. Attractor Dynamics 

If A takes a speciflc value Aq, equation @ integrates 
to give the exponential potential: 



V = Voe- 



-KXocj) 



(5) 



Dark energy with the exponential potential has been con- 
sidered by many authors, e.g. 0, H, The attrac- 
tor dynamics in energy phase-space has been thoroughly 
studied in [28] . Table I summarises the locations and the 
conditions for existence of the attractors. Here we con- 
struct an arbitrary quintessence potential from a series 
of time-slices of exponential potential ([5]) with varying A 
(see [32I, [33I for asymptotic analyses). The dynamical at- 
tractors are thus known exactly at every instant. Hence, 
as explained in detail in the next Section, a large num- 
ber of trial quintessence models may be generated using 
stochastic representation of the function \{z). 



III. GENERATING VIABLE MODELS OF 
QUINTESSENCE 

We shall now impose conditions on the phase-space 
trajectories motivated by present and possible future ob- 
servational constraints. We emphasize that our goal here 
is to gain a physical understanding of quintessence mod- 
els consistent with certain observational constraints. We 
therefore impose representative constraints from future 
experiments rather than modelling the experiments in 
great detail. 



A. Initial Conditions 

We flrst pick an initial point (xo,^o) anywhere in the 
quadrant with uniform probability to start the trajecto- 
ries at high redshift {z 00). We introduce the change 
of variable 



c 



i-C 



(6) 



which maps z G [0, 00) to ^ ^ [O?!)- Thus, initial con- 
ditions at arbitrarily early time can be set at ^ close to 
1. 

We note that an advantage of the phase-space formal- 
ism is that a large (possibly inflnite) range of initial con- 
ditions for quintessence corresponds to a small, flnite re- 
gion in the phase-space. 
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Attractor 


{x,y) coordinates 




Wq 


Condition on A 


Quintessence- dominated 
Scaling solution 


(A/V6, y/1 - AV6) 


1 

3{l + wt)/X^ 


AV3 - 1 

Wt 


A < v^3(l + m) 


(V3/2A2(1 + w,), V3(l - wl)/2X^) 


A > 1/3(1 + Wb) 



TABLE I: Attractors in energy phase-space x—y for the exponential potential ^(0) = Vb exp"^'*'^. The values of the quintessence 
energy density Qq and equation of state Wq at each attractor are also shown. Given the conditions on A in the last column, the 
attractors are either stable nodes or spirals (see for detail). Here, the equation of state of the matter-radiation fluid, Wb, 
is during matter-dominated era, and 1/3 during radiation era. 



B. Interpolating A(^) 

From the initial condition (xo^yo)^ trajectories can be 
evolved at high redshift according to equations (|3|) pro- 
vided the 'roll' variable X{Q is defined at every instant. 
To generate A(C), we assign stochastic values to A in 100 
uniformly distributed C— bins in the range 0.75 < C < 1- 
Values at intermediate ( are obtained by linear interpo- 
lation between bins. 

In each C-bin, we draw A from a uniform distribution 
within the range 

< A < 10. (7) 

Larger values of A correspond to models with steeper 
potentials. Because critical behaviour of the attractors 
only occurs around A = 2 and V^, the range (O is gen- 
eral enough for us to identify a representative behaviour 
of viable quintessence models. A more general scenario 
is where A changes sign, corresponding to potentials that 
are not monotonically decreasing. In Section [VH we ex- 
plain why our conclusions also extend to potentials of 
this type. 



C. Evolution of trajectories 

Trajectories are evolved via Equations ([3|) using the 
C variable until a low redshift {z = 3) where we revert 
to z as the time variable to analyse observables over a 
finer range of bins at low redshift. In z G [0, 3], we make 
a further set of 25 redshift bins between which X{z) is 
interpolated as described above. 

We are interested in the class of models that are 
compatible with tight observational constraints expected 
from future experiments as described in the introduction. 
Let 2;obs denote the redshift at which observational con- 
straints will be tightest. The actual value of 2;obs will 
depend on the range of techniques used to establish the 
constraints. For example, ^^obs ^ 0.3 for supernova sur- 
veys and ~ 0.6 for future weak lensing surveys [2^. For 
this paper, we will adopt ^^obs = 0.3 and note that quali- 
tatively, we can recalibrate to a different value of Zobs by 
applying an additive constant to the redshift. Here, we 
adopt the following constraints: 

(a) the dark energy equation of state Wq is tightly con- 
strained at low redshift to be close to —1 to a per- 



cent accuracy at 2: = Zobs* 

-l<Wq< -0.99 at z = Zohs = 0.3, (8) 
where 



(b) the dark energy density Qq will be tightly con- 
strained at low redshift to a percent accuracy (com- 
patible with the current bounds on Qq from SNIa 
and CMB experiments 0,0): 

0.73 <Qq< 0.74 at z = 0, (10) 

where 

Qq = x^^y\ (11) 

Constraints (a) and (b) define a very small region in 
phase-space through which trajectories must pass z = 
Zohs and terminate at z = 0. Thus, the problem of finding 
viable quintessence potentials translates to finding those 
trajectories that evolve into a small region defined by 
the constraints (a)-(b). This is illustrated in Figure [H 
Trajectories may be straight lines (as in the case of the 
cosmological constant), or complicated paths around the 
quadrant. 

We ran a code that generates a large number of 
quintessence models and selected those that satisfy the 
criteria described in this Section. In this way, we deter- 
mine the generic type of quintessence potentials which 
may be allowed by future dark energy experiments if the 
observations continue to favour w = —1. 



D. Acceptable models 

The trajectories fall into two main categories: 1) static 
and 'skater' models associated with flat potentials, and 
2) dynamical models with trajectories that veer around 
the phase-space in a particular way before terminating. 
We now describe the dynamics of each type of models in 
detail. 



1. Static and skater models 

These models are associated with flat potentials, with 
A :^ throughout their evolutionary histories. Static 
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FIG. 1: Energy phase-space (see Equation ([T])). A viable 
quintessence model can be represented by a trajectory that 
evolve into a small region defined by the imposed conditions 
(Eq. dS} and ^TO\i ). Each point {x,y) corresponds to the ob- 
servables {wq,Qq) via relations Q and (|lip . shown here as 
circular and radial contours. Shaded region denotes where 
cosmic acceleration occurs (it'eff = WqQq < — 1/3). A trajec- 
tory may mimic the cosmological constant by running closely 
along the ^/-axis, or more generally, it may be a complicated 
curved path around the quadrant. 



models have x and behave like a cosmological con- 
stant. For 'skater' models (sj, the quintessence field 
evolves along a flat potential with a constant non-zero ki- 
netic energy. Hence, the phase-space trajectory of skater 
models is the line y = kx, where k = V2V /(p. 

Starting from a generic initial condition within the 
quadrant and given that A <C 1, the trajectories 
quickly evolve towards the quintessence-dominated at- 
tr actor (x^y) ^ (0, 1). To produce static or skater models 
that satisfy our criteria (a) and (b), we found it necessary 
to set the initial coordinates to very precise values which 
were found by backward integration. This simply reflects 
the fine-tuning of the amplitude of the potential if dark 
energy is the cosmological constant. For skater models, 
the constraint (a) on Wq translates to an additional tun- 
ing of the kinetic energy such that 

and so they behave almost identically to the static mod- 
els. 



2. Dynamical models 

The main purpose of this paper is to determine what 
types of dynamical quintessence models will satisfy our 
imposed constraints and to assess whether these are ob- 
servationally distinguishable from the static models. Suc- 
cessful dynamical models evolve in one of the following 
ways: 

(i) The trajectory evolves haphazardly throughout its 
history, avoiding the attractors in Table I. 

(ii) The trajectory evolves towards the scaling-solution 
attr actor, which subsequently decays. 

(iii) If A < 1, the trajectory evolves from a quintessence- 
dominated regime at high redshift to satisfy the 
observational constraints at Zohs- 

It is very difficult to produce type (i) trajectories that 
satisfy the imposed constraints (a) and (b). This is be- 
cause the scaling solution is a global attractor . Tra- 
jectories of type (iii) are unacceptable because an ex- 
tended period of quintessence domination at high redshift 
would have disrupted structure formation. Trajectories 
of type (ii) can produce acceptable models of dynamical 
dark energy, which we now discuss in more detail. 

To satisfy our stringent constraints that Wq is close 
to —1 at ^obs7 clearly all solutions must behave like the 
cosmological constant at low redshift. In the type (ii) so- 
lutions, this behaviour is matched to a scaling solution at 
high redshift, leading to interesting dynamical behaviour. 
Since such a scaling solution has Qq = 3{1 -\- w'^) / and 
necessarily A must be small at low redshift, an accept- 
able model requires a rapid transition of A if the uni- 
verse is to be matter-dominated at high redshift. The 
potentials corresponding to these models therefore have 
a sharp transition from a steep to a shallow slope at some 
characteristic field value. 

Figure [2] shows phase-space trajectories and the cor- 
responding forms of A, Qq and Wq for dynamical models 
that satisfy the imposed constraints. The scaling solution 
attracts trajectories from a wide range of initial condi- 
tions, unlike the static solutions. By construction, all 
of these models have Wq ^ —1 dit z ^ 1. However, un- 
like tracker models considered in the literature, in which 
Wq varies slowly with redshift [36], many models shown 
in Figure [2] display very strong evolution in Wq. The 
evolution of Wq for these models is poorly described by 
parametrizations such diS w = -\- Waz{l z)~^ ^ often 
used in the literature, e.g. [1, [23] and used by the dark 
energy task force to define a figure of merit for comparing 
future dark energy surveys [22]. Note that Wq for these 
models varies wildly above z ^ 5. At z < 3, all of the 
models show very similar evolution in Qq though, as with 
Wq, Qq becomes unconstrained at high redshift. 

To summarise, we have investigated the implications if 
future observations continue to tighten the constraint on 
Wq around —1 with increasingly high precision. Even if 
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z z 

FIG. 2: Examples of trajectories of quintessence models that satisfy the imposed constraints {|8| and (p^ (top panel). From a 
wide range of initial conditions, they approach the scaling solution at early time. But once the scaling solution decays, they 
subsequently evolve towards the region defined by Eq. (|8))- (p^ . As explained in the text, these trajectories are generated by 
stochastic representation of X{z), shown in the second panel (for clarity we have shown only 3 such representations in this 
panel). The lower panels show the variations of the dark energy density Q.q{z) and its equation of state Wq{z), both of which 
show strong evolution at high redshift. 
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FIG. 3: Evolution of Wq and Qq with redshift for the double exponential potential of Equation (O. The parameters have been 
set to = 7 X 10~^^ Vpi, 1^2 = 4 X 10"^^Vpi, A2 = 0.05 with Ai = 25 in panel (a) and Ai = 4 in panel (b). 



future observations were to achieve 1% accuracy in Wq at 
low redshift, our results show that two classes of models 
would be compatible with the constraints: static (skater) 
models that are indistinguishable from a cosmological 
constant, and a class of dynamical models for which Wq 
makes a transition from a scaling solution to Wq = — 1 
at recent epochs. No matter how tight the constraint on 
Wq at low redshift, it is always possible to find dynami- 
cal dark energy models of the second class, which shows 
strong evolution at high redshift. This is the main result 
of this paper. 

IV. RELATION TO THE DOUBLE 
EXPONENTIAL POTENTIAL 

Some of the results of the previous Section can be un- 
derstood by considering the double exponential potential 

V((^) = Fie-^^'"^ + 1/26"^^'"^. (13) 

As discussed above, a pure exponential potential has a 
late time attractor with Wq = Qq = 3{1 -\- Wb) / if 
A^ > 3(1 + Wb) and an attractor with Wq = — 1 + A^/3, 
ftq = 1 if A^ < 3{1 -\- Wb). Evidently, it is not possi- 
ble, using a single exponential potential, to construct a 
model with a matter /radiation dominated phase at high 
redshift with a transition to a quintessence dominated 
phase at low redshift. It is, however, easy to construct 
such models with a double exponential potential ([13]) 
if Ai ^ A2. Figure [3] shows two examples. In Figure 
[3ti, we set Ai = 25, A2 = 0.05, Vi = 7 x lO-^^^m^^, 
V2 = 4: X 10~^^^mpj. These parameters easily satisfy the 
constraint Wq < —0.99 at low redshift, and the high value 
of Ai ensures that the universe is matter /radiation dom- 
inated at high redshift. The potential used to generate 
Figure [3)3 has identical values for Vi , V2 and A2 , but has 
Ai = 4. The lower value of Ai leads to a much stronger 
evolution of Wq with redshift but fails to match the target 



value ofwq < —0.99 by the present day by a wide margin. 
Although the double exponential potential shows some of 
the features of the models described in the previous Sec- 
tion, it is not flexible enough to produce Wq < —0.99 at 
low redshift and rapid evolution of Wq. The numerical 
technique described, however, allows us to find accept- 
able models with both fast and slow evolution of Wq. 



V. DISTINGUISHING BETWEEN STATIC AND 
DYNAMICAL MODELS 

The quintessence potential V{(l)) can be recovered from 
the energy variables by first integrating Equation ([1]) to 
find the quintessence field value as a function of z: 

Hz) = -V6 r ^dz', (14) 
Jo J- + ^ 

where is in unit of reduced Planck mass mpi and is 
arbitrarily set to be zero at z = 0. Next, we integrate 
equation Q to find V{z): 

Viz) = Vo exp (V6 £ ^^f^d.') , (15) 

where Vq is the value of the effective cosmological con- 
stant today. Combining ([T4j) and ([T5]) gives the poten- 
tials V{(j)) as shown in Figure [4] for the models plotted in 
Figure [2l 

The key feature of these potentials is the sudden flat- 
tening at field values characteristic of the transition to 
quintessence domination at low redshift. These models 
can only be distinguished from the cosmological constant 
by observational tests that are sensitive to the steep part 
of the potential at high redshift. Unfortunately, the num- 
ber of feasible observational tests is limited, amounting 
to estimates of distance (either angular diameter or lumi- 
nosity distance) or to measurements of the growth rates 
of fluctuations. 
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and so typically 10^ or more super novae in each of several 
redshift bins would be required to constrain the mod- 
els of Figure [51 Achieving better than a percent accu- 
racy in distance measurements via baryon oscillation sur- 
veys would be formidably difficult even if systematic er- 
rors associated with nonlinearities and biasing were well- 
understood [38]. It seems unlikely, therefore, that dis- 
tance measurements could strongly constrain the models 
of figure [5l 



FIG. 4: The potential y as a function of |A(/)|, the absolute 
change in the quintessence field value with = today. These 
potentials correspond to those quintessence trajectories shown 
in figure [2] The steep parts can be pushed to arbitrarily high 
field values. The prospect of distinguishing these models from 
the cosmological constant is discussed in the text. 



A. Distance measurements 

Luminosity distance dL is given as a function of red- 
shift by 



B. Growth rates 

Growth of large-scale structures can be used to con- 
strain the dark energy equation of state [s^, [10, [U, [42[ . 
Let D denote the amplitude of the growing mode of den- 
sity perturbations. The growth suppression parameter 
f{z) is defined as 



/ 



d\nD 
din a 



The redshift-evolution of / is given by [i^ 



(21) 



dLiz) = —{l + z)x 

-no 



(16) 



dz' 



W{z) = 3 r ^dz\ 
Jo 1 + ^' 



(17) 



The angular diameter distance dA is related to dL by 

dL 



dA 



(1 + ^)2 



(18) 



Figure [5] shows the luminosity distance as a function of 
redshift for the models shown in Figure [H The panel to 
the right shows the fractional change in the luminosity 
distance with respect to a constant A model: 



A^L _ dL{z) - dL,ACDM{z) 



dr 



dL 



ACDM 



(19) 



at low redshift on an expanded scale. One can see that 
over this redshift range, the typical fractional changes 
are of order half a percent. Such small variations in dis- 
tance would be difficult to constrain by either baryon 
oscillation surveys or large surveys of high-redshift type- 
la super novae. For example, if we ignore any systematic 
effects, the fractional error in distance from a sample of 
N supernovae at redshift z with intrinsic magnitude dis- 
persion ai ~ 0.1 is 



dz 



1 



(22) 



where / is normalized to be well into the radiation- 
dominated era. D is then given by 



D{z) = exp 



(23) 



with the normalization D{0) = 1. Figure [6] shows the 
evolution of D{z) for the models of figures [4] and [3 The 
right-hand panel shows the fractional error with respect 
to the growth rate of the constant A model: 



AD ^ D{z) - Dacdm{z) 
D Dacdm{z) 



(24) 



on an expanded scale for low redshift. The variation in 
growth rates is of order a percent or so at low redshift, 
and would be difficult to constrain via observations of 
the mass function of rich clusters 13911 or the integrated 
Sachs- Wolfe effect in the CMB 0,13 • 

In some of the models, the growth rates vary quite 
substantially from the constant A model at high red- 
shift. As discussed in Section IIII D[ this can occur if 
the quintessence potential has a gentle gradient, allow- 
ing solutions in which the dark energy makes a significant 
contribution to the total energy density at high redshift. 
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FIG. 5: Evolution of the dimensionless luminosity distance [Hq / c)dL{z) in the redshift range to 1000 for quintessence models 
generated by the A parametrization as shown in Figured Dashed line shows the corresponding values for the ACDM model. 
The fractional change with respect to the ACDM model (Equation ()19|) ) is shown in the right-hand panel. 



CD 

o 



CD 




FIG. 6: Evolution of the magnitude of growing mode of density perturbation D(z) in the redshift range to 1000 for quintessence 
models shown in Figure 2] with normalization D(0) = 1. Dashed line shows the corresponding values for the ACDM model. 
The fractional change with respect to the ACDM model (Equation ()24p) is shown in the right-hand panel. 



Constraints on this type of behaviour arise from the con- 
sistency of primordial nucleosynthesis (limiting the con- 
tribution of dark energy so that the appropriate expan- 
sion rate at the nucleosynthesis epoch is maintained [45]). 
Additional constraints on the dark energy contribution 
at around the recombination era can be derived from the 
primary CMB anisotropics [8, 46], and in particular, from 
the consistency of the amplitude of the fluctuations at low 
redshift e.g. [47]. Nevertheless, it is clear from Figure [6] 
that it is easy to construct acceptable dynamical models 
in which the dark energy is always subdominant at high 
redshift. 



VI. DISCUSSION AND CONCLUSION 

A. Non-monotonic potentials 

The space of models opens up if we consider 
quintessence in which the field may also roll uphill. In 
this case, trajectories can evolve into the x < quad- 
rant of the energy phase-space. Consequently, trajecto- 
ries that satisfy the imposed constraints ([8|)- (p!Q|) simply 
alternate between the scaling solution in the x > quad- 
rant, and that in the x < quadrant whenever the roll 
direction changes. This gives rise to models which behave 
like the cosmological constant even at high redshift, with 
a fluctuating X{z) of high amplitude. Distinguishing be- 
tween such models and the cosmological constant is even 
more difficult than the monotonic potentials discussed 
above. Thus, we have not considered non-monotonic po- 
tentials in detail here. For specific examples, see (isl, [iot. 
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B. Polynomial potentials and the 'flow' approach 

A number of authors have investigated observational 
constraints on simple low-order polynomial (or Pade ap- 
proximate) parametrizations of the quintessence poten- 
tial [i, lU • In a recent paper [25| , the inflationary flow 
approach of [HO, [HI was adapted to study quintessence, 
but this formulation is equivalent to adopting a low-order 
polynomial for the quintessence potential. 

Such speciflc parametrizations, while useful, cannot 
easily produce the potentials with sharp kinks shown in 
flgure m To illustrate this point, we supplemented the 
differential equations ([3|) by a set of 'quintessence flow 
equations': 



^ = ^VQx\\T^ - 1), 

dTn 

dz 



(25) 



-^VexXFn (1 - Ti - r„ + r„+i) , n > i 

J- ~r Z 



which can be truncated at an arbitrarily high order. 
Here: 



(26) 



is the generalised curvature parameter, with Fi reducing 
to the usual deflnition of F = VV" /V^ usually found 
in analysis of tracker potentials. With a low value of n, 
these equations cannot produce models with a sharp kink 
in the potential. On the other hand, for higher values 
of n, the initial conditions and F^ parameters need to 
be flnely tuned to generate models similar to those in 
Figure [H In contrast, it is straightforward to construct 
dynamical models that satisfy our imposed constraints 
using the A parametrization discussed in Section Hill 



C. Implications 

We have presented a study of quintessence by consid- 
ering stochastic trajectories in energy phase-space (Fig- 
ure [T]). This formalism allows us to generate positive, 
monotonic potentials that are consistent with imposed 
constraints on the dark energy density and its equation 
of state. 



As an illustration of the approach, we have investi- 
gated the implications for dynamical models of dark en- 
ergy should future ambitious observational projects suc- 
ceed in constraining to be —1 at low redshift to an 
accuracy of a percent. 

One set of models that can satisfy these constraints 
is obviously the cosmological constant or 'skater' models 
with extremely flat potentials and a slowly-moving fleld. 
However, we also flnd a second class of acceptable model 
in which the potentials have a sharp transition from a 
steep to a shallow slope as shown in flgure [H For this 
class, the fleld follows a scaling solution at high redshift, 
and so the dynamics is insensitive to the initial condi- 
tions. Although Wq for this type of model can evolve 
strongly with redshift , the deviations in distance mea- 
sures and linear growth rates compared to a constant A 
model over the redshift range < z < 3 are so small that 
it will be extremely difficult to rule them out observa- 
tionally 

Evidently, it is impossible ever to exclude non-trivial 
dynamical models of dark energy no matter how pre- 
cise the observational constraints at low redshift. If, as 
assumed here, the observational constraints on w con- 
tinue to tighten up around —1 to high accuracy, then 
in a Bayesian-Information-Criterion sense, a cosmologi- 
cal constant will be favoured over dark energy models 
involving more parameters [9]. Nevertheless, the numer- 
ical scheme described here shows that it is straightfor- 
ward to construct dynamical models which can satisfy 
observational constraints to an arbitrary precision at low 
redshift, and yet show interesting dynamical behaviour 
at high redshift. Our only handle on this type of model 
comes from constraints on the dark energy density at 
high redshift (for example, from limits on the growth 
rates of ff uctuations [1, \^ ) . Even if the heroic efforts 
of the next generation of dark energy surveys succeed in 
constraining w = —1 to exquisite accuracy, we must re- 
main alert to the possibility that dark energy may have 
been dynamically important at high redshift. 
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